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ABSTRACT 



The time dependent evolution of the accretion disk around black hole is com- 
puted. The classical description of the a-viscosity is adopted so the evolution 
is driven by the instability operating in the innermost radiation-pressure domi- 
nated part of the accretion disk. We assume that the optically thick disk always 
extends down to the marginally stable orbit so it is never evacuated completely. 
We include the effect of the advection, coronal dissipation and vertical outflow. 
We show that the presence of the corona and/or the outflow reduce the ampli- 
tude of the outburst. If only about half of the energy is dissipated in the disk 
(with the other half dissipated in the corona and carried away by the outflow) 
the outburst amplitude and duration are consistent with observations of the mi- 
croquasar GRS 1915+105. Viscous evolution explains in a natural way the lack 
of direct transitions from the state C to the state B in color-color diagram of this 
source. Further reduction of the fraction of energy dissipated in the optically 
thick disk switches off the outbursts which may explain why they are not seen in 
all high accretion rate sources being in the Very High State. 

Subject headings: accretion, accretion disks - binaries: close - instabilities - 
stars: individual (GRS 1915+105) - X-rays: variability 

1. Introduction 

Accreting black holes are strongly vari- 
able sources in all ranges of timescales. 
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Short time-scale variability is usually de- 
scribed through the Fourier analysis and 
the state of the source can be identified by 
a characteristic shape of the power spec- 
tra of galactic black holes (GBH). Physi- 
cal interpretation of this variability is not 
known although some hypotheses were for- 
mulated. Since the power spectra of X-ray 
emission for sources in the Low (hard) and 
High (soft) state are clearly different the 
underlying mechanism may also be differ- 
ent in both cases. Furthermore, the Very 
High State is also characterized by dis- 
tinctively different variability and spectral 
properties so yet another mechanism may 
be responsible for variations in that spec- 
tral state. The nature of this particular 
variability is the subject of our paper. 

Belloni et al. (2000) performed a de- 
tailed observational analysis of the vari- 
ability patterns in the microquasar GRS 
1915+105. It was previously found (Bel- 
loni et al. 1997a, 1997b), that during 
the high luminosity states the spectrum 
of this source is dominated by a ther- 
mal, disk-like component and the variabil- 
ity was modeled phenomenologically by pe- 
riodic disappearance of the inner accretion 
disk (Feroci et al. 1999). Nayakshin, Rap- 
paport & Melia (2000) and Janiuk, Cz- 
erny & Siemiginowska (2000) suggested re- 
cently, that the variability of the micro- 
quasar GRS 1915+105 in its Very High 
State displays features characteristic for 
the limit cycle behaviour well known in CV 
but with the timescales orders of magni- 
tude shorter. 

The limit cycle behaviour of CV is 
caused by the disk thermal and viscous 
instability in the region of partial ioniza- 
tion of the hydrogen (Smak 1984, Meyer & 



Meyer-Hofmeister 1984) and was studied 
in the context of time evolution of these 
objects (Cannizzo 1993, Hameury et al. 
1998). The same instability is also re- 
sponsible for the eruptions of many X-ray 
transients, including GBH, and its char- 
acteristic timescales are days/years (Can- 
nizzo et al. 1995). As the accretion disks 
around supermassive black holes are ex- 
pected to have similar properties, also for 
these objects the analysis of the thermal 
and viscous instability was performed (Lin 
& Shields 1986, Mineshige & Shields 1990). 
The evolution of the accretion disk in AGN 
driven by the partial ionization of the 
hydrogen was computed e.g. in Siemigi- 
nowska et al. (1996). 

However, the standard model of accre- 
tion disk by Shakura & Sunyaev (1973) 
predicts also the existence of additional 
instability due to the radiation pressure, 
with timescales orders of magnitude shorter, 
which is therefore a good candidate to ex- 
plain the behaviour of the microquasar, 
as well as the short term variability in 
AGNs. This radiation pressure instabil- 
ity was first noticed in Pringle, Rees &Pa- 
cholczyk (1973) and studied in Lightman 
& Eardley (1974) and Shakura & Sunyaev 
(1976). 

In order to obtain a limit cycle be- 
haviour in the range of radiation pressure 
instability, some stabilizing mechanism 
must be taken into account. Abramow- 
icz et al. (1988) found that radial advec- 
tion has stabilizing effect on the disk at 
high accretion rates. This was adopted by 
Honma et al. (1991) in the disk evolution 
computations. Thermal stability of the 
transonic accretion disks was studied in 
Szuszkiewicz & Miller (1997) and the limit 
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cycle behaviour was confirmed by calcula- 
tions of Szuszkiewicz & Miller (1998). 

Detailed analysis of the stability prop- 
erties of the disk models in comparison 
with the observational data leads to the 
conclusion that in order to obtain bet- 
ter qualitative agreement with the obser- 
vations, certain modifications of the theo- 
retical models are necessary. Szuszkiewicz 
(1990) studied slim disk models with vis- 
cosity prescriptions different from standard 
a viscosity of Shakura & Sunyaev (1973) 
in order to obtain appropriate outburst 
amplitudes. Also other modified viscos- 
ity laws were introduced (e.g. Milsom et 
al. 1994, Taam et al. 1997). Moreover, 
several authors postulated that ~ 90% of 
the energy generated in the disk is radi- 
ated away by the disk corona (Nayakshin 
et al. 2000) or that energy carried by jet 
has to be included in the model (Janiuk, 
Czerny & Siemiginowska 2000). However, 
while the last modification seems to be well 
supported by the observational data, the 
dominating role of the hot corona is clearly 
in contradiction with the observed spec- 
tral shape of the sources in the Very High 
State. 

In the present paper we propose the 
time dependent model of an accretion disk 
and we calculate the disk evolution based 
on radiation pressure instability. We study 
in detail the model properties and its sensi- 
tivity to the adopted parameters. We con- 
sider the role of advection in the energy 
budget and we take into account the pos- 
sible effects of an outflow and a corona, 
which are consistent with our current un- 
derstanding of the accretion flow at high 
luminosity to the Eddington luminosity ra- 
tio. Finally, we compare the resulting 



lightcurves and spectra with the observa- 
tions. 

2. Radiation pressure instability 

In this section we analyze the local in- 
stability of an accretion disk caused by the 
radiation pressure domination. The sim- 
plest approach to the study of the disk 
stability is based on the vertically av- 
eraged description of the disk structure 
(Shakura & Sunyaev 1973, Shakura & Sun- 
yaev 1976). The vertically averaged model 
offers an easy qualitative insight into the 
disk properties. However, this approach 
is neither accurate nor unique. Various 
approaches to the averaging were used in 
the literature and they differ in detailed 
predictions of the disk stability. 

Paczyhski & Bisnovatyi-Kogan (1981) 
introduced the politropic approximation 
to the disk vertical structure which re- 
sulted in the presence of additional nu- 
merical factors, ranging from 0.5 to 6.0, 
in the averaged equations. The same val- 
ues were later adopted by Abramowicz et 
al. (1988). Honma et al. (1991) adopted 
different coefficients ranging from 1.0 to 
16.0 and in Nayakshin, Rappaport & Melia 
(2000) the coefficients were in the range 
from 0.5 to 2.0. The stability of the disk 
with advection was also studied within the 
frame of the vertically averaged quantities 
(Abramowicz et al. 1988, Szuszkiewicz 
1990, Szuszkiewicz, Malkan & Abramow- 
icz 1996). Here we stress that the values 
of the accretion rate corresponding to the 
transition from the stable branches to the 
intermediate unstable branch may well be 
affected by this approach. 

Recent simulations of the local structure 
of radiation pressure dominated disks per- 
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formed by Agol et al. (2001), Turner & 
Stone (2001) and Turner, Stone & Sano 
(2002) agreed with the analytical models. 
They confirmed that the radiation pressure 
dominated disks are unstable, although de- 
tails of the viscous heating and convection 
may affect the accretion rate at which the 
instability occurs. Interestingly, these first 
simulations also show that the result of the 
instability could lead in some cases to the 
disk collapsing to a new gas pressure dom- 
inated equilibrium. The long term global 
evolution effects of the radiation pressure 
dominated disks have not been calculated 
so far in the 2D simulations. 

Since the full 2D approach to the study 
of global long timescale disk evolution is 
very complex, we adopt the 1D+1D ap- 
proach. Therefore in subsection 2.1 we 
study the full vertical structure of the ac- 
cretion disk with standard prescription 
for the viscosity and advection. Due to 
this approach in subsection 2.2 we are 
able to determine uniquely the appropri- 
ate numerical coefficients, which are sub- 
sequently used in the vertically averaged 
disk model. In subsection 2.3 we include 
in the model the presence of a hot corona 
above the disk and a possibility of a verti- 
cal outflow powered by accretion. In sub- 
section 2.4 we describe the shape of the 
stability curve resulting from the averaged 
equations and we discuss the role of this 
stability curve and of the model parame- 
ters in the time evolution of the accretion 
disk. 

2.1. Vertical structure of a station- 
ary accretion disk 

The disk model at a given radius r is pa- 
rameterized by the mass of the black hole 



M, accretion rate M, and the viscosity pa- 
rameter a. In this work we use a stan- 
dard viscosity prescription, based on the 
assumption that the viscous stress tensor 
is proportional to the total pressure P: 



T T[p = -aP. 



(1) 



While presenting the results, we fre- 
quently use the dimensionless accretion 
rate m, given by 



M 



m 



Mi 



(2) 



Edd 



where M Edd is the critical (Eddington) ac- 
cretion rate: 



M Edd = 



L Ed d AirGMrriH 



c 2 r] 



a T ci] 



(3) 



The efficiency of accretion is r\ = 1/16, as 
it results from the pseudo-Newtonian ap- 
proximation to the disk accretion (Paczyhski 
& Wiita 1980). 

The total energy flux F to t dissipated 
within the disk at a radius r is determined 
by the global parameters: 



3GMM r . , 



(4) 



where f(r) represents the boundary condi- 
tion at the marginally stable orbit: 



f(r) 



1 - 



x3/2 



r — r„ 



(5) 



in the pseudo-Newtonian approximation. 

The disk vertical structure is described 
by a set of basic equations: the equation 
of viscous energy dissipation, 



dF dfl K 



(6) 
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the hydrostatic equilibrium, 
ldP 



p dz 



(7) 



where the total pressure is the sum of the 
gas and radiation pressure: 

P = Pgas + -Prad, (8) 

and equation of energy transfer, 
dT 3kp 



dz 4acT 3 



Fi. 



(9) 



The latter takes into account the presence 
of convection which carries non-negligible 
fraction of energy (see e.g. Shakura, Sun- 
yaev & Zilitinkevich 1978, Dorer et al. 
1996): 

F\ = -^rad V r ad < V a d 

+ F conv V ra d > VadflO) 

Here Fi is the energy flux transported lo- 
cally in the direction perpendicular to the 
equatorial plane, carried either by radia- 
tion or by convection. It differs from the 
flux F due to the presence of radial ad- 
vection. We introduce the advection term 
into our set of equations assuming that the 
fraction of energy radially advected does 
not depend on the distance from the equa- 
torial plane. Such an approach seems to be 
justified, since the dependence of the radial 
velocity on the distance from the equato- 
rial plane is uncertain and it does not influ- 
ence other equations. Therefore we calcu- 
late the required entropy gradients at the 
equatorial plane, denoting the appropriate 
thermodynamical quantities with an index 
e. The local energy flux is then given by 



and the advected fraction / a d v is deter- 
mined from the global ratio of the total 
advected flux to the total viscously gener- 
ated flux (e.g. Paczyhski & Bisnovatyi- 
Kogan 1981, Muchotrzeb & Paczyhski 
1982, Abramowicz et al. 1988): 



/adv 



adv 



tot 



2rF e q SLdv 
'3 Pe GMf(rY 



(12) 



where 



to = (12-10.5A)-^-(4-3ft)-^. 

(13) 

Here (3 is the gas pressure to the total pres- 
sure ratio 



p 



(14) 



Fl — F tot {l — /ady) 



(11) 



Defining the advection term through the 
quantities at the disk equatorial plane is 
an oversimplification which reflects the fact 
that the approach to the vertical structure 
of the disk outlined above does not include 
any information about the vertical distri- 
bution of the radial velocity. The actual 
flow pattern can be studied only if we al- 
low for the departure from the strict verti- 
cal hydrostatic equilibrium. The computed 
flow pattern in Keplerian disks is rather 
complex since the meridional circulation is 
superimposed on the mean accretion flow 
(e.g. Siemiginowska 1988). Therefore, 
there is some uncertainty involved at this 
stage even in the case of disk vertical struc- 
ture included. 

The frequency-averaged opacity k (Rosse- 
land mean) used in vertical structure com- 
putations includes the electron scatter- 
ing as well as the required free-free and 
bound-free transitions. The tables are 
from Alexander, Johnson & Rypma (1983) 
for log T < 3.8, from Seaton et al. (1994) 
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for log T > 4.0 and the value of opacity is 
interpolated between these two tables for 
intermediate values of the temperature, as 
in Rozahska et al. (1999). The numeri- 
cal code which computes the disk vertical 
structure was developed from the version 
from Pojmahski (1986) and subsequently 
modified by Rozahska et al. (1999). 

If the value of the / a d v is specified some- 
how in advance the integration of the disk 
vertical structure proceeds in a usual way. 
The boundary conditions to the set of 
differential equations are specified in the 
equatorial plane z = and at the disk 
surface located at a priori unknown height 
z — H: 

F(0) = 0, F(H) = F tot , 
aT\H)=0.5F tot (l-f adv ). (15) 

The presence of radial derivatives means 
that the problem in general should be 
solved global one, after defining the 
appropriate outer and inner radial bound- 
ary conditions. However, here for simplic- 
ity we adopt the local approach to advec- 
tion developed and successfully used by 
Chen (1995). The numerical value of g a d v , 
given by Equation 13, may be specified at 
each radius by iterations at the two neigh- 
bouring radii. This advection term does 
not depend much on the accretion rate 
and may well be assumed constant on the 
thermal equilibrium curve for a given ra- 
dius. Such a method is correct for disks 
which are Keplerian and not too close to 
the marginally stable orbit, where the con- 
ditions close to the sonic point are essential 
(see also Yuan 1999). 



2.2. Vertically averaged stationary 
disk model 

In the vertically averaged disk structure 
the differential equations (6) - (9) may be 
replaced by the quantities defined as mean 
values or values determined at the equa- 
torial plane. The numerical coefficients C\ 
are uniquely determined here from the ver- 
tical structure model: 



Io H Pdz 
P C H 

J_ ._ Sq P Fdz 
C*2 p e F Ta ,dH 

Jo pzdz 
p e H 2 



Ci = 



(16) 
(17) 
(18) 



Therefore we can use the vertically av- 
eraged relations: 



Ftot = C x -olPSI k H (19) 

is the local energy flux generated through 
the a. viscosity, 



F\H = C*2 



acT 4 
Sup 



(20) 



is the local emitted flux, and the hydro- 
static equilibrium is expressed as: 



C 3 Q 2 K H 2 



(21) 



From the disk vertical structure calcu- 
lations we obtained numerical coefficients 
Ci, Ci and C3 approximately equal to 0.4, 
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2.0 and 0.17 1 (they only weakly depend 
on the accretion rate and viscosity). Us- 
ing these coefficients we can write for any 
given radius the vertically averaged, time- 
independent, energy balance equation: 

(22) 

where the left hand side refers to the vis- 
cous energy dissipation and the right hand 
side describes energy losses due to diffusion 
of radiation and advection. The entropy 
gradient is equal to: 

r^dS P , , 

T—r = —q adv (23) 
dr p 

where g a d v is described by the Equation 
(13) and is approximately constant. The 
energy balance Equation (22) together 
with the averaged hydrostatic equilibrium 
Equation (21) allow us to compute the sta- 
bility curve on the surface density versus 
temperature (X — T) plane. These sta- 
tionary solutions are presented in Section 
2.4. 

2.3. Stationary disk with a corona 
and vertical outflow 

Detection of the power law tail in the 
hard X-ray spectra of black hole X-ray 
binaries and AGNs (Mushotzky, Done & 
Pounds 1993, Tanaka & Lewin 1995) in- 
dicates, that somewhat hotter medium is 
present in the vicinity of a standard accre- 
tion disk. Here we assume, that this hot 

Our factors d correspond to the coefficients de- 
fined in Muchotrzeb & Paczyhski (1982): B x = 
0.8, B 3 — 5 and B 4 — 6, while these authors use 
B x = 0.5, B 3 = B A = 6 and Honma et al. (1991) 
use Bx = 1.0, B 3 = 16 and B A = 8. 



gas forms a corona above the disk, which 
extends down to the marginally stable or- 
bit and dissipates a constant fraction / cor 
of the total accretion energy. This energy 
does not contribute to the disk black body 
emission, but is stored in the corona and 
radiated in the hard X-rays. 

For a number of objects the presence of 
a collimated vertical outflow is confirmed 
(e.g. Mirabel & Rodriguez 1994 for GRS 
1915+105). We also include in our model 
such a possibility, however we do not adopt 
any collimation mechanism. The plasma 
ejection is powered by the energy released 
in the innermost disk parts. Therefore 
we assume that another fraction of a to- 
tal energy is converted into the wind and 
is not seen in the disk multicolor black 
body emission. We parameterize this frac- 
tion using the luminosity to the Edding- 
ton luminosity ratio, so that for accretion 
near the Eddington limit we should expect 
the outflow to be the strongest. Unlike 
in the ADAF models, where a powerful 
wind may take away a considerable amount 
of mass and angular momentum (Bland- 
ford & Begelman 1999), we assume, that 
in our disk the carried mass and angular 
momentum may be negligible (see Discus- 
sion). Therefore the presence of a wind is 
reflected only in the energy equation. 

The energy balance Equation (22) is 
modified by the outflow and the corona. In 
principle, in the time independent model, 
we could put the outflowing and coronal 
fractions either in the cooling or in the 
heating term with no change in the final 
energy balance. However, in the physical 
case both the outflow and the corona play 
their role in cooling of the accretion flow, 
so in general the emitted flux will have the 
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form: 



1 nrT 4 

F- = C 2 -^—[(l - / cor )(l - f out )]-\ 



where 



H 3np 



/out — 1 — 



1 



1 + Am 2 



(24) 
(25) 



we take after Nayakshin et al. (2000). Our 
model parameters are A and f cor , which 
we assume to be constants, ranging from 
to 1. Below we describe the role of these 
parameters in the resulting stationary disk 
solutions. 

2.4. The local stability curve 

When the energy balance given by 
Equation 22 is satisfied, the stationary lo- 
cal disk solutions form a characteristic S 
curve on the E — T plane (alternatively 
on the vertical axis we may plot accretion 
rate M instead of the temperature). In 
this section we consider the curves calcu- 
lated for the 1OM black hole. In the case 
of a supermassive black hole the disk tem- 
perature is ~ 2 orders of magnitude lower 
and the surface densities may reach ~ 10 
times higher values for a given Rj Rs c hw It 
results in shifting the entire stability curve 
down and right on the E — T diagram. 

The accretion disk is thermally and vis- 
cously unstable when ^ < (Lightman & 
Eardley 1974, Meyer & Meyer-Hofmeister 
1981). This is the case for the accretion 
rates in the range (mi,m 2 ), where the S 
curve has a negative slope. The lower 
branch of the curve represents stable solu- 
tions, resulting from the gas pressure dom- 
ination. The intermediate unstable branch 



appears due to radiation pressure domina- 
tion, while the upper hot branch is sta- 
bilized by advection (Abramowicz et al. 
1988). The upper and lower transition 
points rhi and rh 2 depend on the global 
model parameters (i.e. M and a) as well 
as on the local distance r. The depen- 
dence on the global parameters is relatively 
weak, but there is a strong dependence on 
the radial location r: a change of r results 
mostly in a strong shift of the entire curve 
up or down on the diagram. It translates 
into strong dependence of the extension of 
the radiation-pressure dominated disk re- 
gion on the accretion rate (e.g. Hure 1998, 
Janiuk & Czerny 2000). 

In Figure 1 we plot the solutions for 
several disk radii. If the external accre- 
tion rate does not exceed a critical value, 
the whole disk remains in the cold, sta- 
ble regime and the matter is transported 
inwards with constant velocity. However, 
for higher external accretion rates the in- 
ner disk parts may find themselves in an 
unstable mode. Then the disk undergoes 
oscillations, heating up and cooling down 
in a cycle. This may cause the observed lu- 
minosity outbursts or flickering. The range 
of accretion rates for which the disk is lo- 
cally unstable varies across the disk. It is 
wider at small radii and becomes narrower 
at larger radii. This effect is even stronger, 
if the full vertical structure is computed 
instead of the constant C\ coefficients in 
the vertically averaged model. In this case 
at large radii, r > 200-R Sc hw the whole in- 
termediate branch of the S curve gets a 
positive slope and the instability vanishes. 
As a result the outer parts of the disk are 
stable, while the inner parts may undergo 
oscillations. 
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Fig. 1. — The stability curves calculated 
for M = 10M Q , and 4 different radii : R = 
fOi?schw (solid curve), R = 60i?schw (dot- 
ted curve), R = 140i? S chw (short-dashed 
curve) and R = 290i?s chw (long-dashed 
curve). The viscosity parameter a = fCT 2 . 
The effects of outflow and corona are ne- 
glected. Note, that on the vertical axis 
we plot local dimensionless accretion rate, 
m = M I ' M-Edd, instead of the disk temper- 
ature. 



The critical accretion rate, above which 
the disk becomes unstable, as well as the 
outbursts amplitude depend on the loca- 
tion of the turning points m\ and m 2 on 
the stability curve. This is determined by 
the presence and the strength of the disk 
hot corona and/or outflow. 

In Figure 2 we plot the stability curves 
at location R = I0i?s c hw, resulting from 
different assumptions about the fraction 
of energy dissipated in the corona. When 
the corona is neglected, the lower turning 
point corresponds to a very small accretion 
rate (rhi ~ 0.06) and the outburst ampli- 
tudes should be very large. The presence of 
the corona increases the critical accretion 
rate and reduces the size of the unstable 
branch, which results in a reduction of the 
outbursts amplitude. 

In Figure 3 we plot the solutions ob- 
tained for different strengths of the out- 
flow. The strong outflow suppresses the in- 
stability at high temperatures and reduces 
the amplitude of an outburst. However, it 
does not affect the cold stable branch and 
has no influence on the lower turning point 
on the stability curve. 

The evolution time scale and duration 
of an outburst is determined by the vis- 
cosity. Here we use the constant viscos- 
ity parameter a. In Figure 4 we plot the 
dependence of the disk solutions on this 
parameter. The stability curve is shifted 
almost horizontally with a. Therefore the 
viscosity parameter does not influence out- 
burst amplitudes, while its value may make 
them last shorter or longer. 
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Fig. 2. — The stability curves calculated 
for M = 1OM , R = lO.Rschw and 4 values 
of the fraction of energy dissipated in the 
corona: f COT = 0.0 (solid curve), f cor = 0.3 
(dotted curve), / cor = 0.5 (short-dashed 
curve) and / cor = 0.9 (long-dashed curve). 
The viscosity parameter a = 10~ 2 . The 
outflow effect is neglected. 



Fig. 3. — The stability curves calculated 
for M = 1OM , R = 10i? S chw and 4 val- 
ues of the outflow parameter A (see Equa- 
tion 25): 0.0 (solid curve), 0.025 (dot- 
ted curve), 0.125 (short-dashed curve) and 
0.325 (long-dashed curve). The viscosity 
parameter a = 10~ 2 . The effect of the 
corona is neglected. 
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3. Time dependent disk 
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3.5 



log(otE) 



Fig. 4. — The stability curves calculated 
for M = 1OM , R = 10R Schw and 3 
values of viscosity parameter: a = 1CT 1 
(solid curve), a = 1CT 2 (dotted curve), and 
a = 1CT 3 (dashed curve). The corona and 
outflow effects are neglected. 



In this section we describe the time de- 
pendent disk model. Here we treat the 
disk radial structure independently from 
the vertical one, which is included in our 
calculations via the coefficients C\ in the 
averaged equations. The crucial param- 
eter, which links the radial and vertical 
disk structures, is the viscosity parameter 
a. The viscosity mechanism is responsi- 
ble for both radial transfer of the material 
inwards and for the local dissipation of en- 
ergy, which in turn is either radiated verti- 
cally or advected. Once a is specified, we 
adopt the vertically averaged equations in 
the time evolution of the disk. 

3.1. The basic equations 

The energy balance expressed by the 
Equation (22) gives the stationary solution 
in the E — T plane. In order to construct 
the time dependent model we solve the two 
basic equations for S(r, t) and T(r,t), fol- 
lowing the evolution of the accretion disk. 
The mass conservation equation: 



~dt 



1 dM 
27rr dr 



where 



M = — 27rrEv r , 



(26) 



(27) 



is combined with the angular momentum 
conservation equation: 

M^(r 2 fi) = - J^rX)- (28) 

Here Q is the Keplerian velocity and the 
np component of the stress tensor is equal 
to T rip = aPH = (3/2)fti/£. P is the total 
pressure, H is the disk height, v denotes 
the kinematic viscosity and v T is the radial 
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velocity. From these equations we obtain 
the equation of the surface density evolu- 
tion: 

<9S 1 d / , „ d 



dt r dr 



(3rV^ (r V Vs) ). (29) 



The disk temperature evolves due to the 
departure from the thermal energy conser- 
vation: 



Erf = «-«" 



(30) 



Here we adopt viscous heating 

Q+ = C 1 (3/2)QaPH (31) 
and radiative cooling 

Q- = C 2 [(l-/ cor )(l-/out)]- 1 (4 ( 7/3/ t )(T 4 /S). 

(32) 

The coefficients C\ and C 2 were defined 
in Section 2.2. Combining Equation (30) 
with time-dependent forms of Equations 
(23) and (13) we obtain the energy equa- 
tion: 

<91nT <91nT 



+ v r - 



dt dr 
4-3/3 /91ns d\nH 



12 - 10.5/5 V & 



dt 



+ 



5lnS\ Q + -Q~ 

+ V <—) + (12-10.WPH - (33) 
The radial velocity is given by: 

v r = "f ^JWn (34) 

The corona and the vertical outflow, pa- 
rameterized by / C or and / out (see Section 
2.3), play a crucial role in our calculations. 
In principle, they could contribute either 
to the cooling or to the heating of the disk, 



and therefore influence strongly the results 
of the time evolution. Here we choose 
again the physical case, e.g. corona and 
outflow as the disk cooling mechanisms, 
consistently with our stationary solution. 
In the following we assume that they add 
to the cooling Q~ , while the heating term 
Q + includes only viscosity. The corona pa- 
rameter f COT is assumed constant with time 
and radius, while / out depends on time and 
radius of the disk via the local accretion 
rate rh(r,t); in this case the constant pa- 
rameter of the evolution is A (see Equation 
25). 

3.2. Numerical method 

We solve the equations (29) and (33) via 
explicit method, using the constant time 
step, comparable to the thermal time scale 
(Wallinder 1991): 



.r 1 / 2 M 

itherm = 8.75 X 1(T 5 (r - 1)— [s] 

a M G 

(35) 

where r is the radius expressed in Schwarzschild 
radii. 

The viscous time scale depends on the 
disk thickness: 



^visc ~ 



S 1 /r 

dt 



^therm 

^therm 
tvi^r 



f ^ 

r / 



(36) 



For typical values of H/r in our model the 
viscous time scales are 3-4 orders of mag- 
nitude longer than the thermal time scales. 

The equation (29) is transformed into a 
simple diffusion equation due to the conve- 
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nient change of variables (Smak 1984): 

dZ _ 12 d 2 
~dt ~ y*dy* 
y = Zr 1 ' 2 E = yY, (37) 

On the uniform radial grid in y we com- 
pute in every time step the new values of 
the 6 variables: 2™, T™, P™, Hp, (5f and 
up, where j is the number of radial zone 
and n is the number of a time step. Initial 
values are taken from the stationary solu- 
tion on the lower branch of the S— T curve. 
The boundary condition on the inner disk 
radius, 3P S chw, is that T in = S in = 0. The 
outer disk radius is chosen to have a stable 
solution, determined by the external, con- 
stant accretion rate m cxt This assumption 
implies, that a certain amount of the angu- 
lar momentum is taken away at the outer 
radius of the computed disk region. The 
typical radial extension in our calculations 
was 300i?schw and the number of points in 
the radial grid was 220. 

3.3. Results 

In this Section we present the results 
of the time dependent calculations. The 
main parameters of the model are mass of 
the central black hole M, external accre- 
tion rate m ext , viscosity a, coronal fraction 
/ C or and outflow strength A. Throughout 
this work we assumed the mass of the cen- 
tral object equal to 10 M , which is ap- 
propriate order of magnitude for galactic 
black hole binary systems. The viscosity 
was constant in the disk and in every model 
we adopted a = 0.01. 

The external accretion rate is the ac- 
cretion rate outside the region of the ra- 
diation pressure instability, at the outer 



radius, which in most cases was fixed at 
300 -Rschw This accretion rate value was 
chosen in the range m cxt = 0.19 — 0.71, 
so that the inner disk parts were unstable. 
However, we did not take into account any 
effects connected with possible ionization 
or gravitational instability, which may be 
present further out in the disk. We per- 
formed our calculations for the full range 
of corona and outflow strengths (/ cor and 
A ranging from to 1). 

Once all the model parameters are spec- 
ified, the disk starts to evolve, and tries to 
reach the state specified by m ext . At every 
disk radius the solution is found along the 
lower branch of the stability curve S, up to 
the moment when the innermost radii find 
themselves in the unstable mode. Then the 
oscillations begin, and the disk tempera- 
ture and density change periodically with 
time. 

Below we present the results of our cal- 
culations. Firstly, we discuss the evolution 
at the particular disk radius during one cy- 
cle, e.g. starting and ending on the lower 
stable branch of the S curve. This trans- 
lates into one outburst of the disk total 
luminosity. Secondly, we show, how the 
luminosity of the disk changes with time, 
depending on the assumed corona and out- 
flow efficiencies. Finally, we present the in- 
fluence of the radiation pressure instability 
on the geometrical thickness of the accre- 
tion disk, and we show an example of the 
disk spectral energy distribution. 

3.3.1. Evolution of the surface density 
and temperature during one cycle 

In Figure 5 we show the evolution of the 
surface density during one cycle. Between 
the outbursts the radial dependence of E 
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is smooth and has a positive slope. The 
radial inflow through the accretion disk is 
stable and the material accumulates slowly 
in the inner region. When the instabil- 
ity starts, the density drops down rapidly 
in the innermost radii, along with the in- 
creasing local accretion rate. The mat- 
ter is transferred inwards much faster than 
it is supplied from the outer disk regions 
and the inner disk empties. At the outer 
end of the unstable region there appears a 
local density peak, which shifts outwards 
as the instability propagates into larger 
radii. The outburst finishes when there is 
no more additional matter supply. As the 
disk cools down, the density in the inner 
part gradually rises and when the cycle is 
completed, the flow returns to the initial 
configuration. 

In the outburst the disk temperature 
rises mostly in the inner disk regions. In 
Figure 6 we show the temperature evo- 
lution during one cycle. The instability 
begins at the innermost disk radius. It 
propagates outwards and after several time 
steps the heat front reaches further re- 
gions. The temporal decrease of the tem- 
perature in the inner part is characteristic 
for models with both high external accre- 
tion rates and considerable corona share. 
During this temporal decrease the disk is 
still heated up in outer regions, and this 
phase corresponds to the smaller "dip" in 
the lightcurve. It is followed by the sub- 
sequent heating of the entire unstable disk 
region, which is the cause of the sharp lu- 
minosity increase. At the end of the cycle 
the whole disk cools down rapidly. The 
cooling front again starts from the inner 
disk radius and propagates on the thermal 
time scale. The disk enters into the phase 



of relative faintness of the inner parts. It 
may observationally create an impression 
that the inner disk has been completely 
evacuated, while actually the disk always 
extends down to the marginally stable or- 
bit. 

The evolution on the £ — T plane de- 
pends on the local disk radius. The inner- 
most parts of the disk reach the highest 
temperatures and if the external accretion 
rate is high enough, they may even jump to 
the upper, advection supported branch of 
the stability curve. This depends on other 
model parameters, e.g. the strength of the 
corona and the outflow. Further out the 
upper branch is never reached and the disk 
undergoes oscillations in the intermediate 
unstable region. The outermost parts of 
the disk always remain on the lower, cold 
stable branch (see Figure 7). 

3.3.2. The luminosity of the disk 

The total disk luminosity is equal to 
L = j R ° ut 2nrF TSu ddr and it results only 
from the radiation of the disk. We do 
not include the fraction radiated by the 
corona, carried by the outflow or advected. 
The coronal fraction is assumed constant 
during the evolution. However, the out- 
flow strength depends on the luminosity to 
the Eddington luminosity ratio and varies: 
from negligible value between the out- 
bursts to / out ~ 0.35 during the outburst. 
Therefore, in the absence of an outflow, 
the resulting outburst amplitude remains 
constant, being somewhat lower for a disk 
with a strong corona than for a disk alone. 
When the outflow is present, the observed 
outburst amplitude varies during the evo- 
lution. The observed luminosity fraction 
radiated by the accretion disk is equal to 
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fobs = (1 - /cor) (1 - /out)- Depending on 
other model parameters, in the case of a 
strong outflow its luminosity may reach 
even / out / / obs = 30 - 70%. 

In Figures 8 and 9 we show the lumi- 
nosity curves for m cxt = 0.7. In the first 
case the outflow was not included and the 
fraction of the energy dissipated in the 
corona was equal to / cor = 0.5. This 
resulted in strong and regular outbursts. 
For the model with the outflow contribu- 
tion parameterized by A = 0.075 the out- 
burst shape is quite irregular. In this case 
the share of the corona has been dropped 
to / C or = 0.3, which allowed us to keep 
the mean outburst amplitude roughly un- 
changed. 

In order to check if the luminosity vari- 
ations would vanish for the largest outflow 
and corona shares, we computed a grid 
of models for the full range of these two 
parameters. Because the outflow creates 
irregular outburst shapes and amplitudes 
we parameterized the luminosity variations 
using the root mean square approximation, 
defined as: 



rms 



<L> 



where 



2 j t ; i {L{t)-<L>fdt 
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(38) 
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(40) 



The time range (t\, £2) was chosen arbitrar- 
ily and in our case it is (1.4, 1.41) x 10 5 
sec, long enough to cover several outbursts 
for each of the studied lightcurves. In 
Figure 10 we plot the luminosity varia- 
tions in terms of rms for different corona 



and outflow strengths. Both outflow and 
corona suppress the variations. Therefore 
the largest outbursts amplitude is achieved 
for disks with no corona and a very weak 
outflow, while the strong outflow and/or 
corona results only in luminosity "flicker- 
ing" , which could be hardly detectable. 

3.3.3. The disk thickness 

The geometrically thin disk approxima- 
tion is always satisfied in our model. How- 
ever, during the outburst the inner parts 
of the disk become thicker than in the cold 
state (see Figure 11). This is because the 
temperature achieved during the outburs 
is controlled by the upper, advection sup- 
ported branch of the stability curve. 

3.3.4- Radiation spectrum 

In Figure 12 we show the radiation spec- 
tra emitted by the accretion disk in the 
maximum (outburst) and minimum lumi- 
nosity phases. We have made a simplifying 
assumption, that the disk radiates locally 
as a black body of temperature T e s(r) and 
the resulting spectrum is integrated over 
the entire disk. We adopt the temperature 
distribution in the unstable region as that 
computed from our time dependent model, 
and we assume that the outside disk is sta- 
ble and has a standard temperature distri- 
bution T c fr ~ r - 3 / 4 . xhe spectrum between 
the outbursts is a standard disk black- 
body spectrum. However, when the disk 
becomes unstable, the radial temperature 
distribution is changed, which affects the 
spectral slope. Also, the maximum tem- 
perature is much higher, which causes the 
shift of the spectrum towards higher fre- 
quencies (see Zampieri et al. 2001). The 
presence of the corona and/or the out- 
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flow should add a power law component 
to the total spectrum, with an extension 
and slope appropriate for the hot plasma 
properties. This spectral component is not 
plotted in the figure. 

4. Discussion 

4.1. General properties of the model 

The radiation pressure instability oper- 
ating in the innermost parts of the accre- 
tion disk may be responsible for the peri- 
odical luminosity changes of a number of 
sources. The limit cycle behaviour pro- 
duces outbursts with timescales and am- 
plitudes which could reproduce the ob- 
served variability of Galactic black hols (on 
timescales of minutes) and short term vari- 
ability of AGNs (on timescales of years), 
provided that certain assumptions about 
the model are satisfied. 

In our time dependent calculations we 
use vertically averaged equations of the 
disk structure. Appropriate numerical co- 
efficients need to be used in the vertical 
averaging, because they influence the crit- 
ical accretion rate for which the limit cycle 
starts working. The description of the disk 
vertical structure also influences the com- 
puted outburst amplitudes, particularly in 
the case of the 'disk + corona' model. 

We calculated the time- dependent disk 
model and associated lightcurves and spec- 
tra assuming that a moderately strong 
corona extends above an accretion disk 
and contributes to the total accretion en- 
ergy. Our assumption is consistent with 
the observed hard X-ray tail in the spectra 
and soft to hard X-ray luminosity ratios 
(e.g. Zycki, Done & Smith 2001 for soft X- 
ray transients, Janiuk, Czerny & Madejski 



2001 for the quasar PG121 1+143). We also 
take into account the presence of the verti- 
cal outflow. These two assumptions let us 
obtain good agreement with the observed 
GBH lightcurves. For a bare accretion disk 
without a corona or an outflow, the out- 
bursts are much too strong in respect to 
observations and the maximum luminosity 
exceeds the Eddington limit. 

The outflow efficiency suppresses the 
maximum luminosity during the outburst. 
The luminosity reaches even 0.504 L Edd 
in the case of no outflow (model with 
/cor = 0.5), while for the outflow strength 
A = 0.01 the luminosity drops to 0.2LEdd 
(for the same corona share). The total out- 
flow power could in principle be as high as 
that emitted in the disk (Falcke & Bier- 
mann, 1999). Fender (2001) shows that the 
outflow is likely to be comparable in power 
to the integrated X-ray luminosity. In our 
calculations the outflow power is variable 
and depends on L/L-^^ ratio, achieving 30 
- 70% of the observed disk luminosity dur- 
ing the outburst, for various sets of param- 
eters. 

Not only the outflow but also the corona 
reduces the outburst amplitude. If the 
fraction of energy dissipated in corona is 
small, a strong outflow is needed to sup- 
press outbursts: in case of / cor = 0.1 
the outflow efficiency parameter has to be 
A > 0.25. For f cor = 0.5 only a small 
amplitude flickering (5logL < 0.05) occurs 
when A > 0.05. Here we adopt a simpli- 
fied description of a corona, assuming con- 
stant coronal dissipation fraction f COT , but 
in general this parameter may depend on 
the radius (Janiuk & Czerny 2000). Re- 
cently, Merloni & Fabian (2002) proposed 
that the coronal fraction may drop for high 
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accretion rates, being the most efficient in 
the low state of the source. 

When either an outflow or a corona are 
strong enough (/ cor > 0.55 and/ or A > 
0.25), the outbursts do not occur and the 
disk produces only flickering. Therefore 
these two mechanisms support each other 
and the stronger the corona the weaker 
outflow is required to suppress disk's out- 
bursts. We illustrate this correlation in 
Figure 13, where we show on the A — f cor 
plane the parameter range for which we de- 
tect significant outbursts. Here the crite- 
rion for outburst detection was chosen that 
rms > 0.03, that is the mean luminosity 
variations are above ~ 5% in logarithm. 

The results strongly depend on the 
adopted external accretion rate. For small 
m ex t the whole disk remains stable and the 
outbursts do not occur. However above a 
certain critical value the instability devel- 
ops in the innermost disk part. The critical 
accretion rate above which outbursts are 
expected is m cxt = 0.28 for a disk without 
a corona, but it increases to m cxt = 0.68 
when the corona shares 50% of the en- 
ergy. This critical accretion rate does not 
depend on whether an outflow is present 
or not (see Figure 3). The amplitude of 
the outburst only weakly depends on the 
external accretion rate. Instead, the dura- 
tion of the bursts increases strongly with 
m ext , and the separation between the out- 
bursts increases as well. For instance, for 
fhext = 0.56 the burst duration is St ~ 250 
seconds, while for rh ext = 0.39 it is only 
St ~ 150 seconds (model with no corona). 
It is the result of the enlarged instability 
region when the accretion rate is high. 

The outburst amplitude is determined 
by the separation between m\ and m 2 , e.g. 



the lower and the upper turning points on 
the S curve (see Section 2.4). The point 
rhi is controlled mainly by whether the 
corona is present or not and how large is 
its contribution to the total accretion en- 
ergy. The point rfi2 is influenced mainly 
by the strength of the outflow. In order to 
keep the outburst amplitudes no too high 
and not exceed the Eddington limit in the 
maximum luminosity phase, the point m 2 
must never be too high. In principle, the 
small value of wi2 could alternatively be 
achieved by changing the viscosity param- 
eter a, and would not require any contri- 
bution from the hot wind. The variable 
viscosity may be parameterized e.g. by 
the disk thickness: a ~ (h/r)~ n . How- 
ever, we found that for low viscosities (a < 
10~ 3 ) the unstable radiation pressure dom- 
inated branch becomes too small for the 
instability to develop, and finally the disk 
achieves quasi-stable state in the interme- 
diate range of temperatures (see Wallinder 
1991). 

In our model we assume the a vis- 
cosity prescription of Shakura & Sunyaev 
(1973). If the magneto-rotational instabil- 
ity was taken into account, the a param- 
eter would also contain a magnetic contri- 
bution, which could slightly increase its ef- 
fective value (Armitage et al. 2001). The 
qualitative results however should not be 
changed in this case. On the other hand, 
we assume a stationary material supply 
to the radiation pressure instability zone 
from the accretion disk regions outside ~ 
200-Rschw! while in principle additional in- 
stability mechanisms, such an ionization 
or gravitational instability, may play cru- 
cial role further out in the disk (Menou & 
Quataert 2001 and references therein). 
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In all computed cases the duration of 
the outburst was significantly shorter (by 
a factor of 3 or more) than expected on 
the basis of the viscous timescale computed 
at the outer edge of the instability strip 
from the analytical formulae provided by 
Shakura & Sunyaev (1973). Therefore, nu- 
merical computations are needed to give 
appropriate viscous time scale estimates. 

4.2. Application to GRS 1915+105 

Observations of the microquasar GRS 
1915+105 outbursts imply rather short 
time scales and low amplitudes of an out- 
burst. Ferocietal. (2001) show lightcurves 
of the microquasar obtained by Beppo 
SAX. Typical soft X-ray band outbursts 
with AlogL ~ 0.6 amplitude last for 
At ~ 100 sec. The disk luminosity in 
the outburst never exceeds Eddington lu- 
minosity. Therefore, based on our model, 
the whole limit cycle should run between 
rfii ~ 0.05 and wi2 ~ 1.0. 

The values of the transition accretion 
rates rhi and m 2 (see Sections 2.4 and 
4.1) depend significantly on the adopted 
description of the disk structure. In par- 
ticular, the value rrt\ resulting from the 
model of Nayakshin et al. (2000) will be 
higher by a factor of 5 when the correct 
factors computed from the vertical struc- 
ture model are included. It weakens signifi- 
cantly the motivation to postulate a strong 
coronal dissipation (up to 90% in Nayak- 
shin et al. 2000) in this source, which is not 
observed but which served to increase rhi 
by a factor of 10. Therefore in our model 
we propose a moderately strong corona, 
dissipating 30-50% of the total energy. If 
the outflow takes away required amount of 
the total emitted power depending on the 



L/L-Edd ratio, the maximum observed disk 
luminosity should be lowered. 

In our calculations we adopt the black 
hole mass equal to 10M Q (recent near in- 
frared spectroscopy results indicate the 
mass of 14 ± 4M , Greiner 2001) and 
constant viscosity parameter a = 0.01. 
With a corona dissipating 50% of the to- 
tal energy, no outflow, and external accre- 
tion rate equal to 70% of Eddington rate 
the model results in strong and regular 
outbursts of AlogL = 0.7 amplitude and 
timescales of At ~ 140 sec. If the outflow 
is present than the required by the ob- 
servations corona is weaker and dissipates 
only 30% of energy. In this case the the- 
oretical lightcurve is quite irregular, with 
shorter and longer outbursts lasting for 
At ~ 50 — 100 sec and variable amplitudes 
of AlogL ~ 0.3 - 0.7. 

The scenario of radiation pressure in- 
stability has one advantage over the disk 
evacuation scenario (Belloni et al. 1997b). 
As stressed by Rao, Yadav & Paul (2000), 
the rise/fall timescale during outbursts, 
t ~ 10s, is too short for a cold disk to be 
removed or reinstalled, as this should hap- 
pen on a viscous timescale. In our model 
the whole accretion disk exists all the time, 
just its thermal state changes during out- 
burst, with thermal timescale being short 
even for relatively high accretion rates. We 
see that clearly from Figure 2 of Janiuk, 
Czerny & Siemiginowska (2000) - rise and 
fall are much shorter than the duration of 
long outbursts typical for high accretion 
rates. 

Belloni et al. (2000) pointed out, that 
the source may occupy three distinct re- 
gions on the color - color diagram in X- 
rays. The colors are defined as HRi = 
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F(5 - l?,keV)/F{2 - hkeV) and HR 2 = 
F(13-60keV)/F(2-5keV). The state A 
corresponds to small values of both HR\ 
and HR 2 , state B is for HR\ > 1.1 and 
small -ffi?2, while state C is for small HRi 
and large i?i?2- The source alterations be- 
tween states A, B and C are observed pe- 
riodically, however the only one way tran- 
sition between states B and C has been 
observed, namely from B to C, while the 
reverse transition has never happened. In 
our model we could identify the state C 
with the disk occupying the low stable 
branch of the S curve (see Figure 14), with 
a very soft disk blackbody emission and 
dominating power law component from the 
hot corona. The state A corresponds to 
the turning point on the stability curve, 
when the contribution of the disk black- 
body emission is substantial, and state B 
refers to the unstable branch on the S 
curve. In the last case the emission from 
the hot disk is more pronounced in the X- 
ray spectrum. The disk normally evolves 
along the stability curve, going from C to 
A and further to B, or back to C state, de- 
pending on the accretion rate. However, 
when the disk reaches the unstable region 
B, the only possibility is to drop down to 
the state C, when the system finds the 
shortest possible way to come back to the 
stable solution. 

4.3. Luminosity states in other sources 
and their relation to radiation 
pressure instability 

Galactic black holes (GBH) exhibit a 
number of distinct luminosity states (van 
der Klis 1994, Poutanen, Krolik & Ryde 
1997). Low (hard) state is characterized 
by a power law spectrum, with the photon 



index r ~ 2, and a complex variability. 
High (soft) state spectra are dominated 
by the disk emission, and the power spec- 
trum is a simple power law with the index 
~ 1 and a high frequency cut-off above 100 
Hz. When the source gets brighter, in the 
Very High State, the spectrum is complex, 
with considerable and frequently variable 
contribution of the power law and disk 
emission, and the power spectrum shape is 
closer to the one in the hard state. Finally, 
the faintest sources represent a quiescent 
state. 

Our model based on the radiation pres- 
sure instability predicts that the luminos- 
ity threshold exists and the disk becomes 
unstable when the luminosity exceeds this 
threshold. There is an important question 
to be put forward: are the Low, High and 
Very High states related in any way to the 
radiation pressure instability? 

The theoretical luminosity threshold de- 
pends on the corona efficiency and varies 
from 0.28L Edd for no corona (/ cor = 0) to 
0.68L Edd for f COT = 0.5. The range of ob- 
served accretion rate corresponding to var- 
ious luminosity states is rather difficult to 
assess. Firstly, only the luminosity in a re- 
stricted energy band is usually available, 
while estimates of the bolometric luminos- 
ity are reliable only if the data truly repre- 
sent the broad spectrum. Secondly, there 
may sometimes be a confusion in the state 
definitions. The Very High State looks 
similar to both the hard state and the soft 
state, and it may be classified as one of 
those two, unless a careful analysis of both 
the spectral and temporal behaviour is per- 
formed. Keeping this in mind, we have 
collected a number of accretion rate esti- 
mates from the literature for the purpose 
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of qualitative discussion. In Table 1 we 
show the estimates of the accretion rate 
(or luminosity to the Eddington luminosity 
ratio). Information is not complete since 
many states are observed mostly in tran- 
sient sources and then the evolution may 
be too fast to collect many data (or some- 
times even a single data set) for a given 
luminosity state. 

There seem to be a large difference be- 
tween the range of accretion rates in differ- 
ent states among various objects. System- 
atic errors and dispersion in inclination an- 
gles may partially be responsible for those 
observed differences, but there may also be 
some additional parameters hidden in the 
description of formation of the hot (coro- 
nal?) plasma. Still, sources are in their 
Low States when rh ~ 0.05, in the High 
States when m ~ 0.1, and in the Very 
High States when m ~ 0.5. The luminos- 
ity range corresponding to the High state 
is quite narrow, factor of 2 - 3 in accretion 
rate between the minimum and maximum 
value. 

The Low and the High states are ba- 
sically expected to be stable. The differ- 
ence between these states does not seem 
to be related to radiation pressure instabil- 
ity but there must be a mechanism of hot 
plasma generation which works less effi- 
ciently when the total flux is larger. There 
are several phenomenological models pro- 
posed but the physics of formation and ge- 
ometry of hot plasma is under discussion 
(e.g. Poutanen 2001, Di Salvo et al. 2001). 
Some objects in the High State are possi- 
bly slightly above the threshold predicted 
by our model although their coronal emis- 
sion is not stronger than / cor = 0.5 (LMC 
X-3 and J1748-288). It should be consid- 



ered whether these sources are indeed in 
the High State, or rather in a disk domi- 
nated Very High State. 

Sources in their Very High State are 
definitively in the range corresponding to 
the radiation pressure instability. Out- 
bursts lasting a hundred to a thousand 
of seconds are expected. Such outbursts 
are however seen only in microquasars - 
GRS 1915+105 (see Sect. 4.2) and GRO 
J1655-40 (Remillardet al. 1999). The lack 
of outbursts in other sources is consistent 
with our model only if the corona and/or 
outflow are strong enough to suppress the 
instability completely. Since the 'corona' 
(observationally - a power law component) 
in the Very High State contributes typi- 
cally to about 30 % of energy, considerable 
outflow is also needed (parameter A in our 
model of order of 0.2). In sources, which 
do not exhibit jets, the outflow must be 
uncollimated. The model does not include 
any description of the jet formation pro- 
cess, therefore it basically applies to such 
uncollimated outflow. 

It remains to be checked whether this 
outflow is consistent with observations. 
Energetically, the request for such an out- 
flow is plausible since the observed colli- 
mated jet may carry a considerable frac- 
tion of a total energy (Fender et al. 1999, 
Reeves et al. 2001). However, this mate- 
rial, if uncollimated, is in our line of sight 
and covers the disk, so part of the 'corona' 
effect in the Very High State may be due 
to the Comptonization in this outflowing 
material. Such winds were also discussed 
in several contexts (e.g. wind driven insta- 
bility of Shields et al. 1986, Beloborodov 
1999, ADIOS models of Blandford and 
Begelman 1999). In the present model 
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the wind should develop as a result of the 
disk instability, thus providing a stabilizing 
mechanism, and therefore it should cover 
only the radiatively unstable part of the 
disk. This hot wind plasma should have 
different properties than the hot coronal 
plasma formed in the Low state, which is 
probably produced by another mechanism. 

4.4. The vertical outflow 

In this work we postulate that a certain 
amount of the accretion energy is taken 
away by the vertical outflow. Since we 
do not include any collimation mechanism, 
this outflow has a form of a hot wind rather 
than a collimated jet, however its origin 
are the innermost regions of the accretion 
flow. We did not study this outflow in de- 
tail, since we took into account only the 
energy loss due to the wind. We neglected 
the mass loss and the effect of the hot flow 
onto the disk spectrum. However, we can 
estimate the properties of the outflowing 
wind a posteriori. The energy flux of the 
outflowing material is given by: 



be given by: 



out 



(41) 



where v z is the outflow velocity in the ver- 
tical direction and m z is the mass loss rate. 
Therefore the outflow efficiency, defined by 
the Equation (25), is equal to: 



out 



—. x — 2 , 42 

3M/(r) V 



where local flux F tot is given by Equation 
(4) and is the Keplerian velocity on the 
orbit r. Therefore the mass loss rate will 



m 7 



w (1 + ^) x ^» vw < 43 > 

and may be negligible if e.g. v z > Sv^,. 

For the exemplary model with A = 0.1, 
we set rh ~ 1 and (f v /f z ) 2 ~ 0.1 at 
r = 10-Rschw- The estimated density of 
the wind is then equal to p w j n d = iii z /v z = 
2.7 x 10~ 9 g cm -3 . This means that the op- 
tical depth is of the order of t cs ~ 0.03 and 
the outflowing medium should not have 
any substantial impact on the observed 
spectral properties of the X-ray source. 
The Comptonization in the wind may lead 
at most to the formation of a weak, steep 
power-law tail, additional to the disk black 
body emission. The X-ray reflection ampli- 
tude is also reduced by the bulk motion in 
the range of /^buik ~ 0.6 — 0.7, so this effect 
is practically negligible (Janiuk, Czerny & 
Zycki 2000). 

We stress that the uncollimated out- 
flow which partially stabilizes the radia- 
tion pressure dominated disk is launched 
by a different mechanism than the jet seen 
in GRS 1915+105 during the prolonged C 
states. The latter is observed in the low 
luminosity intervals and is probably con- 
nected with the radio emission (Klein- Wolt 
et al. 2001). 

5. Conclusions 

Time dependent evolution of the inner 
accretion disk is based on the radiation 
pressure instability, playing major role in 
so called 'a disks'. This process accounts 
for a number of properties of Galactic black 
hole binary systems. 



• The radiation pressure dominated in- 
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Table 1 

Luminosity ranges in Eddington units 



Object 


Very High State 


High State 


Low State 


References 


J1748-288 


0.30 - 0.50 


0.05 - 0.25 


0.01 - 0.04 


Revnivtsev et al. 2000 


Nova Muscae 


0.09 - 0.50 


0.03 - 0.09 


0.006 


Zycki et al. 1998 


LMC X-3 


_b 


0.07 - 0.3 


< 0.07 


Wilms et al. 2001 


Cyg X-l 


_b 


> 0.03 


< 0.03 


Gierlihski et al. 1997, 










Gierlihski et al. 1999 


GS 2023 


< 1.6 


0.07 


0.04 - 0.05 


Zycki et al. 1999 a,b 


GRO J1615 


0.2 - 0.25 


0.1 


a 


Remillard et al. 1999, 










Gierlihski et al. 2001 


LMC X-l 


_b 


0.15 


_b 


Gierlihski et al. 2001 


GX 339-4 


0.5 


a 


0.05 - 0.15 


Miyamoto et al. 1995, 










Wilms et al. 1999 


GRS 1915+105 


> 0.15 


_b 


_b 


Zdziarski et al. 2001, 










Belloni 1997b 










Trudolyubov et al. 1999 


a No observations in this state 
b No transitions to this state 
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nermost regions of the accretion disk 
are unstable and may oscillate be- 
tween cold and hot phase. When 
the hot front goes through the inner 
disk, its surface density drops down 
rapidly by a factor of 3. As the 
disk cools down, it switches to the 
state of a relative faintness and the 
X-ray spectrum may be dominated 
by a power law emission from a hot 
corona. However, the disk is never 
evacuated completely and extends al- 
ways down to the marginally stable 
orbit around a black hole. 

• For high enough accretion rates 
(rh > 0.28) and moderate strengths 
of the outflow and the corona, the 
instability results in a characteristic 
variability pattern, with 100 - 200 
sec. outbursts of an amplitude reach- 
ing AlogL ~ 0.3 — 0.7. The irregu- 
lar outbursts' shape, observed in the 
microquasar GRS 1915+105, may re- 
sult from the presence of the outflow, 
which takes away a high fraction of 
the energy when the disk luminosity 
is close to the Eddington limit. 

• The presence of luminous outbursts 
and their amplitude depends on the 
assumed strengths of the disk corona 
and the outflow. These two elements 
may weaken the outbursts or even 
suppress them completely, therefore 
in most systems accreting at high rh 
(observed in the Very High State) a 
rapid, large amplitude variability is 
not detected. 

• The distinct regions occupied period- 
ically by the microquasar on the X- 
ray color-color diagram correspond 



in our model to the characteristic 
phases during the time evolution. 
The theoretically predicted evolu- 
tionary track through these states: 
C — > A — > B — > C explains the 
lack of C — > B transitions. Observed 
transitions B — > A and A — > C and 
the fast X-ray flickering are not ac- 
commodated yet in the model, since 
they are both probably caused by hot 
plasma instabilities. 
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1 2 

log (RAscw) 

Fig. 5. — The surface density evolution 
during a single outburst. Model parame- 
ters: M = 1OM , a = 10" 2 , m cxt = 0.7. 
The corona share is f COT = 0.5. The outflow 
is neglected. 



5.5 I 

1 2 

log (RAscw) 

Fig. 6. — The temperature evolution dur- 
ing a single outburst. Model parameters: 
M = 1OM , a = 10"2, m ex t = 0.7. The 
corona share is / cor = 0.5. The outflow is 
neglected. 
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log (2) log (2) 



Fig. 7. — The evolution on the £ — 
T plane for four different values of ra- 
dius: 5.75i? S chw, 9.88.Rschw, 30.5i? S chw and 
140.5i?schw The solid line marks the lo- 
cal stability curve and the points represent 
the solutions of the time dependent disk 
model. We show the time sequence from 
the beginning of the evolution through a 
few outbursts. Model parameters: M = 
1OM , m cx t — 0.65, a = 0.01. The corona 
share is neglected and the outflow effi- 
ciency parameter is A = 0.025. 
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Fig. 8. — The disk lightcurve. Model pa- 
rameters: M = 1OM , rriext = 0.7, a = 
0.01. The fraction of the energy dissipated 
in the corona is / cor = 0.5. The outflow is 
neglected. 



28 




100 



200 300 
time [s] 



Fig. 9. — The disk lightcurve. Model pa- 
rameters: M = 1OM , m ex t = 0.7, a = 
0.01. The fraction of the energy dissipated 
in the corona is f COT = 0.3 and the outflow 
parameter is A = 0.075. 




Fig. 10. — Disk luminosity variations 
for different outflow and corona strengths. 
The open triangles mark solutions for 
corona parameter f COT = 0, open squares 
for / C or = 0.1, open circles for / cor = 
0.2, filled triangles for / cor = 0.3, filled 
squares for f COT = 0.4 and filled circles for 
/cor = 0.5. Other model parameters were 
kept constant: M = 10M Q , m ex t = 0.7, 
a = 0.01. 
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Fig. 11. — The H/R ratio during the disk 
evolution in the minimum (solid line) and 
maximum (dashed line) luminosity states. 
Model parameters: M = 1OM , m cxt = 
0.7, a = 10~ 2 . The corona share is f COT = 
0.5. The outflow is neglected. 



Fig. 12. — The disk blackbody spectrum 
in the minimum (solid line) and maximum 
(dashed line) luminosity states. Model pa- 
rameters: M = 1OM , m ext = 0.7, a = 
0.01. The corona share is / cor = 0.5, but 
the power law tail is not plotted here. The 
outflow is neglected. 
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0.1 0.2 0.3 

A 

Fig. 13. — Parameter range, for which the 
outbursts are expected. Model parame- 
ters: M = 1OM , m ext = 0.7, a = 0.01. 
The shaded region marks the corona and 
outflow parameters, for which the luminos- 
ity variations are greater than 5% (rms > 
0.03). 
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Fig. 14. — The hardness ratio states A, B, 
and C, defined as in Belloni et al. (2000), 
identified with the different positions of the 
disk on the stability curve. The curve is 
plotted for 10 -Rschw The cycle performed 
during the time evolution is marked by ar- 
rows. 
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